Landau damping in thin films irradiated by a strong 
laser field 

D F Zaretsky 1,3 , Ph A Korneev 2 *, S V Popruzhenko 2 3 
and W Becker 3 

1 Russian Research Center "Kurchatov Institute", 123182, Moscow, Russia 

2 Department of Theoretical Nuclear Physics, Moscow Engineering Physics 
Institute, Moscow, 115409, Russia 

3 Max-Born-Institut, 12489, Berlin, Germany 

Abstract. The rate of linear collisionless damping (Landau damping) in a 
classical electron gas confined to a heated ionized thin film is calculated. The 
general expression for the imaginary part of the dielectric tensor in terms of the 
parameters of the single-particle self-consistent electron potential is obtained. For 
the case of a deep rectangular well, it is explicitly calculated as a function of the 
electron temperature in the two limiting cases of specular and diffuse reflection 
of the electrons from the boundary of the self-consistent potential. For realistic 
experimental parameters, the contribution of Landau damping to the heating of 
the electron subsystem is estimated. It is shown that for films with a thickness 
below about 100 nm and for moderate laser intensities it may be comparable with 
or even dominate over electron-ion collisions and inner ionization. 



Submitted to: J. Phys. B: At. Mol. Opt. Phys. 
PACS numbers: 



* e-mail: korneev@theor.mephi.ru 



Landau damping in thin films irradiated by a strong laser field 



2 



1. Introduction 

The interaction of strong laser fields with nanobodies such as clusters and thin films has 
become a subject of wide interest over the last decade. A very significant distinguishing 
feature of clusters and thin films irradiated by a strong laser field as compared with 
a gas of isolated atoms and even with solids and plasmas is the very efficient energy 
transfer per atom from the laser field to the respective matter [1] . The high efficiency 
of energy absorption in nanobodies produces fast heating of their electron subsystem 
up to temperatures on the keV scale and leads to the occurrence of many interesting 
new phenomena such as Coulomb explosion of clusters, production of ions with very 
high charges, x-ray emission from cluster targets and films, to name a few. Therefore, 
the mechanisms of energy absorption by the electron subsystem in nanobodies are 
of primary interest and importance for understanding many aspects in strong-field 
nanophysics. 

The physics of laser-irradiated nanostructures including clusters and thin films 
have been reviewed recently El EH ■ Various models of the energy absorption from a 
laser pulse were discussed. In the nanoplasma model |Sj , which is meant to describe 
clusters irradiated by a strong laser field with intensities of I0 15 -I0 18 W/cm 2 , the 
imaginary part of the dielectric constant e = s' + ie" is supposed to be the result 
of individual electron- ion collisions like in an infinite plasma, so that e" — 2v e i/u>. 
Here v c \ is the frequency of electron-ion collisions and w the laser frequency. In a 
plasma with parameters typical for intense-laser-cluster interaction, the time between 
individual collisions is around 10-100 fs. For laser pulses of moderate intensity (10 13 - 
10 14 W/cm 2 ), as they are used in studies of the optical breakdown in dielectrics (see 
El EJ and references therein), electron-phonon collisions give a relaxation time on 
the picosecond scale. 

Another efficient mechanism of energy absorption is vacuum heating 9 . However, 
it is relevant only for a not too thin plasma and for high laser intensities when the 
conditions a > xq > I are satisfied, where a is the thickness of the plasma, xq the 
amplitude of the electron oscillations under the action of the laser field, and I the size 
of the skin layer. For the parameters considered here, namely films with a thickness 
a < 100 nm and moderate laser intensities, the relation between a, xq and I is rather 
the other way around, so that the mechanism of vacuum heating is irrelevant. 

Finally, laser energy may be absorbed in electron collisions with the confining 
surface. This effect can be interpreted as collisionless damping, or Landau damping, 
in a spatially restricted plasma. Such a mechanism of energy transfer in various 
multiparticle systems from collective into individual-particle degrees of freedom was 
first considered for longitudinal waves in a quasi- neutral infinite plasma |1UII11| . For a 
finite system, Landau damping was studied for the first time in 1966 |12j for a collection 
of cold fine metallic particles, assuming an external field weak enough to leave the 
electron gas at zero temperature. In a cold spherical nanoplasma, the imaginary part 
of the dielectric constant produced by Landau damping has the form 

e" = e l (uj) = C(u) — , (1) 

where vf is the Fermi velocity, a the spatial size of the system (e.g. a — R for 
a spherical cluster), and the dimensionless function C(oj) depends on the boundary 
conditions for the electrons and on the geometry of the system. A detailed discussion 
of Landau damping in cold metal clusters including an analysis of the role of boundary 
conditions for the electrons and a comparison with experimental data may be found 
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in references [131 1141 ITS ) . For cold particles with radii below 100 nm, collisionless 
damping dominates over collisions with phonons |14| . viz. £l 2> £oa, so that the 
absorptive part (JIJ determines the cross sections of light scattering and absorption by 
a small nanobody with a cold degenerate electron subsystem. 

For nanostructures interacting with a strong laser field, the state of the electron 
gas differs significantly from that in cold clusters: first, the electron gas is heated to 
a temperature on the keV scale, so that the electron plasma becomes purely classical; 
second, a nanobody irradiated by a strong field becomes partially ionized so that the 
self-consistent electron potential may differ considerably from the rectangular well 
used for the description of cold particles |121 1131 H5) . However, the contribution of 
Landau damping to absorption remains substantial |161 1171 I18| . In particular, it 
defines the amplitude of the electric field inside a cluster near the Mie resonance 
|17| . The latter may become particularly important, when it is nonlinearly excited. 
This was recently proposed theoretically ^3 EDI an d observed experimentally |2J m 
laser-irradiated clusters. In references the width of the Mie resonance was 

estimated assuming that the frequency of electron collisions with the boundary is 

^5' (2) 
where v is the thermal or effective electron velocity. This crude estimate yields for the 
imaginary part of the dielectric constant the expression Q with the effective velocity 
v in place of Upt an d C = 2. Obviously, it does not take into account the shape of 
the self-consistent potential whereby the electrons are trapped nor the shape of the 
thermal distribution function. Therefore, the effect of Landau damping in a heated 
classical nanoplasma is still lacking a quantitative treatment. 

In the present work, we calculate for the first time the contribution to the dielectric 
tensor from collisionless damping in the heated electron plasma trapped inside an 
ionized laser-irradiated one-dimensional (ID) thin film. The plasma is supposed to be 
classical, so that this contribution does not depend on Planck's constant. As we shall 
show, in a classical finite system the simple form |J2j is only realized under specific 
conditions. In general, the behavior of the imaginary part £l of the dielectric tensor 
is quite complicated and depends on the shape of the self-consistent single-particle 
electron potential and on the specific boundary conditions. Its numerical value can 
be substantial, so that this mechanism may be dominant in the process of electron 
plasma heating. Some preliminary results of this work were annonced in |22| . 

The paper is organized as follows. In the next section, we state the problem and 
obtain the general expressions for the rate of collisionless energy absorption and for 
the dielectric tensor e,fe (ui) in a homogeneous isotropic thin film. The case of diffuse 
reflection of electrons from the boundary of the self-consistent potential is considered 
in the third section. The fourth section contains the application of the results obtained 
to the description of laser light absorption in a thin film. Details of the calculations 
are delegated to appendixes A and B. 

2. Statement of the problem and general expression for the dielectric 
tensor in a thin film 

Let us consider a one-dimensional system (a thin film) consisting of ions 
homogeneously distributed inside a layer of width a and a heated electron plasma 

f In 1171 . the effective velocity was taken in the form ti cf f = \Jv^. + v^, where vt and v q are the 
thermal and quiver velocities, respectively. 
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concentrated mainly within this ion layer. As a result of the partial evaporation of the 
electrons, the system has some nonzero positive charge and produces a self-consistent 
static electric potential U(z), which traps the residual electrons. We treat the ionic 
subsystem as frozen, which usually is a valid assumption on the femtosecond time 
scale!. The self-consistent potential obeys the Poisson equation 

MJ{z) = -4ne[ziri?(z) - n° e (z)}, (3) 

where e is the absolute value of the electron charge, nf(z) and n®(z) are the equilibrium 
concentrations of ions and electrons, respectively. Inside the ion layer, i.e. for 
—a/2 < z < a/2, we have nf(z) = n . Equation © is not sufficient to define U(z), 
because the equilibrium electron distribution n®(z) is unknown, too. For the consistent 
definition of these two functions, equation must be solved simultaneously with the 
Boltzmann equation for an electron plasma. Here, we do not consider this problem 
and take the self-consistent potential as a given function. 

Oscillations of the electron cloud in such a system arise either due to the action 
of the external laser field or due to a small deviation of the electron cloud from its 
equilibrium position. We suppose that the rate of electron-electron and electron-ion 
collisions is low, so that oscillations of the electron cloud are far off equilibrium and do 
not change the electron distribution in momentum space. The oscillations produce the 
perturbation Sn e (x,t) of the spatial electron distribution. If their amplitude is small 
with respect to the spatial size a of the layer, the oscillations may be considered as 
homogeneous (dipole) with the laser frequency oj§. This approximation corresponds 
to the model of an incompressible fluid for the electrons. For clusters in strong laser 
fields, this model was used in reference and its applicability was supported by 
numerical simulations |2()j . The time-dependent electron density then is given by a 
rigid displacement by from its equilibrium density n®(z): 

n c {z, t) = n° e (z) + 5n c (z, t) = n° c (z - £(*)), (4) 

where the displacement obeys the nonlinear oscillator equation | 

i + 2 7 HC + W *F(0 = cosM + a), (5) 



2k-i 

(6) 



fc=i 

where the external force is provided by the laser fields 

E(t) = E cos(wi + (7) 

and we define the z axis perpendicular to the plane of the film (see Fig. 1). The 
damping constant 7 in equation J5J is related to the imaginary part of the dielectric 
tensor [see equation <(Tol below]. 

I For clusters irradiated by a strong laser field, the time when this assumption ceases to be valid has 
been estimated, for example, in reference 1231 . This time can be measured experimentally via the 
dependence of the intensity of the scattered light on the laser pulse duration Both the estimates 
and the experimental results support the conclusion that for nanobodies with radii around 100 A the 
approximation of frozen ions holds up to a few hundred femtoseconds. 

§ Free oscillations take place at the plasma frequency u) p = 47re 2 nore/m c , where m c is the electron 
mass and re is the fraction of electrons trapped inside the ion layer. 
| For a spherical cluster, such equations are derived in references 1251 1191 1261 

% Of course, the phase a has no physical relevance at this point. It could be removed by just resetting 
time. However, it will be convenient in the treatment of diffuse reflection below. 
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Here we consider linear damping. Therefore, it is enough to retain the first term in 
the expansion of F in equation 10. The resulting harmonic oscillator has the solution 

£(*) = eE ° z cos(cut + a + /3), 

m e \ (ujp ~ uj 2 ) 2 + 47 2 cj 2 

2 ^ (8) 

J(UJ 2 ~ LU 2 ) 2 + A"f 2 LU 2 

The electric field inside the thin film is the sum of the external field and the 
induced field. Assuming 7 « w we obtain 



E(£, t) = 47reno£e 2 + E cos(ujt + a) = 6q± cos(ujt + a) + e z £ z cos(uit + a — S), 

£ ± = E _l, £ z = = = , sin<5 = ( —) sin ,3.(9) 

^(1-^/^2)2+4^2 \ujJ 

The rate Q of collisionless energy absorption is equal to the work per unit time of 
the electric field (JJJJ on the individual electrons averaged over their distribution. For 
a classical system as we consider it here, it is independent of Planck's constant. Yet, 
surprisingly, the simplest route for its calculation is via quantum mechanics. Following 
24 , we consider electrons populating energy levels in the symmetric potential well 
U{z) — U(—z), with some distribution function p(e) and subject to the uniform electric 
field @. The electron energy e depends only on the principal quantum number n. 
The field produces electron transitions from their initial level n to a final level n\ . The 
rate of energy absorption from the laser field per unit time and per unit square of the 
film is 

Q = hujn c a P(tn)(wfi, ni - w n,nJ' ( 10 ) 

where w~ ni and w+ ni are the transition probabilities per unit time of an electron 
from the level n to the level n\ owing to absorption or emission of a laser photon, 
respectively. The distribution function p may also depend on time, due to evaporation 
of electrons from the well. Moreover, individual collisions or fluctuations of the self- 
consistent potential may influence the energy absorption. However, here we restrict 
ourselves to the simplest limiting case, when individual collisions and fluctuations may 
be neglected and the quantum state of the electrons may be described in terms of wave 
functions. Classically, such an approximation corresponds to the motion of electrons in 
a given time-independent potential with specular reflection from the boundary whose 
position depends on the electron energy. The opposite limiting case of diffuse reflection 
of the electrons from the boundary will be briefly considered in the next section. 

We calculate the probabilities w^~ n in first-order perturbation theory with respect 
to the field 10- Although the external laser field is strong, far from the plasma 
resonance the amplitude of electron cloud oscillations satisfies Xq ~ eEQ Z /m c uj 2 <C a. 
So, for a laser pulse with intensity I ~ 10 15 W/cm 2 , a film with thickness a ~ 10 nm, 
and the electron density n c ~ 10 23 cm -3 , we have oj p s» 8 eV and xq 1 nm. These 
estimates support the applicability of lowest-order perturbation theory. The details 
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of the calculation are given in Apendix A. Finally, then, the rate of energy absorption 
has the form 



2 r-2 00 

nn e ae 

k=0 



|p2fc+i(efc)| 2 d 

2m~^ |n' e (e fc )| de 

Here the sum extends over all roots of the equation 
(2A; + l)n(e fc ) = w, 



mm 



(ii) 



(12) 



where 17(e) is the round-trip frequency of an electron with energy e in the self-consistent 
potential U(z). The condition l|12|) admits a simple and transparent interpretation: 
a nonzero contribution to energy absorption comes from those electrons only whose 
period of motion in the potential U(z) is equal to an odd number of laser periods. 
Therefore, as in the case of an infinite plasma or cold finite plasma, energy absorption 
proceeds through "resonant" particles. 

The expression for the dielectric tensor £ik{ui) may be found by a comparison of 
the result l|ll|) with the standard formula for the rate of energy absorption from a 
plane electromagnetic wave interacting with a dielectric |28j , 



Q = -^ReP -E*, 



(13) 



where Pi = (e^- — 8ik)£k/^ is the polarization per unit volume, £ is the electric field 
(J3J inside the film, and E the incident field J7J. For a homogeneous isotropic film, 
the dielectric tensor should be proportional to the unit tensor 5^. However, Landau 
damping produces anisotropy since only oscillations along the z axis lead to energy 
absorption. As a result, the dielectric tensor in a homogeneous isotropic film has the 
form 

, ,2 - 



£ifcM = £oM + 1 1 )Sik + ie L (^)S lz S k 



(14) 



where £0 = + i^o is the contribution from the polarizability of the bound electrons 
and electron- ion collisions. A nonlinear contribution to e'q due to multiphoton 
ionization can also be included. The anisotropic term £l is responsible for Landau 
damping. 

Expanding the external field into p-polarized and s-polarized waves as shown in 
Fig. 1, we get from equations l|13|l and (|14fl (cf. figure 1 for the definition of s and p) 



n - au} p 2 



"0 J2 



p z sin^ + £o (s +P cos 2 1?) 



Comparing equations (JTTJ, l|15|). and © we find 



7H = 



2u) 

IT 



2 00 



E 

fe=0 



b2fc+i(e fc )| 2 d 
\Q' e (e k )\ de 



(15) 

(16) 
(17) 



Expression 11711 is the main result of the present work. It describes the 
contribution of Landau damping to the dielectric tensor for a homogeneous isotropic 
thin film. 
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3. Dependence on the boundary conditions 

In our calculations so far, we have assumed that an electron moves inside a 
thin film under the action of the field Q and the self-consistent single-particle 
potential U produced by the ion core and the quasi-equilibrium distribution of all 
electrons. Individual collisions and fluctuations both have been completely ignored. 
In macroscopic language, this approximation corresponds to the specular reflection of 
an electron from the boundary of a self-consistent potential. Below, we identify the 
dielectric constant calculated in this approximation by the superscript "s", £l — > ££■ 
In this case, the influence of the electric field JjJJ on the electron motion is a coherent 
effect accumulated over many periods. This leads to the feature expressed by equation 
(JI2J): the energy absorption from the electric field is possible only for selected electrons 
with "resonant" energies as explained below equation l)12fl . All other energy levels 
do not contribute to e£. However, the approximation of specular reflection is only a 
limiting case, since collisions are present and do have an influence on the efficiency of 
energy transfer. The reason is that they interrupt the coherent accumulation of the 
phase of the electron motion. In consequence, energy absorption is possible even if 
the resonance condition (|12|l is not satisfied. This may result in increased absorption 
in comparison with the case of specular reflection. 

A rigorous quantum-mechanical account of collisions requires a density-matrix 
approach. However, an estimate of the effect of collisions on Landau damping may 
be obtained using a simple classical method. The influence of individual collisions 
on the electron trajectory in a self-consistent field increases when the velocity of the 
electron is small, i.e. in the vicinity of the turning points. Here, an electron can 
experience a large momentum transfer in an individual collision and, as a result, will 
lose the memory of its phase. To estimate the effect of this phase mismatch, we use a 
simple model for diffuse reflection of an electron on the boundary. We suppose, that 
between the boundaries every electron moves in a self-consistent potential under the 
action of the field J^J, but in the collision with the boundary the phase is reset, so 
that subsequent trips of the electron between the boundaries are independent. In this 
approximation, the total energy that the electron collects from the field © is given 
by the incoherent sum of the contributions from the various segments of the electron 
motion between the turning points. Treating the electron classically, instead of (f 1 01) 
we have for the rate of energy absorption by the ID system 



where A(e, r) is the energy absorbed during the time r by the electron with the energy 
e averaged over the initial conditions. Obviously, this energy is the incoherent sum 
of the contributions from the various halfperiods within the time r [see equation (|19|l 
below] . 

In the linear approximation used in this paper, averaging over the initial 
conditions is equivalent with averaging over the phase a of the field at the instant 
of time when the electron starts from the boundary, so that 



where r(e, a) is the time the electron with energy e needs to travel across the well from 
one boundary to the other, provided it started its motion at t = when the electric 




(18) 




(19) 



o 



o 
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field has the phase a. This time differs from a half period of the unperturbed motion, 
viz. 7r/f2(e), by terms that are linear in the amplitude of the field. The quantity N is 
the number of half periods within the time r. so that N — r/(7r/0(e)). 
The velocity v(e, t) of the electron is to be found from the equation 

1 dU e - 

v z = E z cos(ut + a) (20) 

m c oz m c 

with initial conditions z(0) = Zo(e), where zo(e) is one of the turning points for the 
motion at the energy e. 

In contrast to the case of specular reflection when the general expression (| 1 1 ft 
for the rate of energy absorption is, in principle, a functional of the self-consistent 
potential, equation (121 III cannot be solved in closed form for an arbitrary potential 
U(z) even in linear approximation with respect to the field (j!|Jl + - Therefore, the 
case of diffuse reflection can be considered analytically only for some particular self- 
consistent potential. In the next section, we will consider and compare the results for 
specular and diffusive reflection for the simplest case of an infinitely deep rectangular 
well. 

4. Discussion 

If the distribution function /(e) of the electrons and the self-consistent potential U(z) 
arc known, the contribution to the dielectric constant Ijl7|) and all relevant quantities 
such as the reflection and absorption coefficients and the rate of energy absorption 
can easily be calculated from equations Q11JI- (|14|) . and i|17fl . either numerically or 
even in analytical form. A separate and still unsolved problem is the determination 
of the quasi-equilibrium distribution function and the self-consistent potential for 
electrons in a heated ionized film. Here we consider the arguably simplest model: the 
infinitely deep rectangular ID well, which permits a complete analytical treatment 
both for specular and for diffuse boundary conditions. This is an example of a 'closed' 
potential, where electrons cannot leave the system even for arbitrarily high energy. In 
other words, there is no evaporation of electrons from the film after the self-consistent 
potential has been established. In such a potential, global thermodynamic equilibrium 
will be achieved and electrons will be distributed according to the Gibbs law, so that 

This model provides a simple analytic expression for the dielectric constant, which 
allows one to understand general properties of collisionless absorption in a classical 
finite system. 

Details of the calculations are given in Appendix B. Both for specular and 
diffuse boundary conditions the dielectric constant £l is the product of the factor 
Wp/w 2 and a function of the single parameter x = KT / irm c uj^ > a 2 = v^/uj v a. This 

feature simplifies the analysis of the behavior of the rate. Plots of the function £l 
calculated numerically from equations 1|B.2(1 . i|B.8() . and l|B.9(l are given in Fig. 2. 
In all cases, the low-temperature behavior of the dielectric constant is linear with 
respect to x ~ vt, i.e. the simple law is recovered with in place of vf and 

+ This means that the solution of equation 121 it cannot be analytically expressed via the solution 
of the unperturbed equation as was done in the previous section for specuilar reflection. In linear 
approximation with respect to the field J^J, equation 1201 reduces to the Hill equation. 
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with the dimensionless function C(oS) — CqU> 2 /lj 2 , where Cq waries between 1 and 
2 depending on the boundary conditions* (see insert in Fig. 2). In the opposite 
case of high temperatures, the behavior of the dielectric constant is different for 
specular and diffuse reflection. However, realistic experimental parameters (see the 
estimates below) usually correspond to the limit i < 1, when the result is practically 
insensitive to the boundary conditions. The behavior of the dielectric constant versus 
temperature may be easily understood from the resonance condition l|12|) . Indeed, 
the oscillation frequency of a particle in the rectangular well increases with increasing 
energy as [see IjB.lll ]. This means that the upper 'resonant' level has the energy 
m e 0J 2 a 2 1 2tt 2 [corresponding to k = in the resonance condition (|12|) ]: all other 
'resonant' energy levels are at lower energies, near the bottom of the well. As a 
result, for low temperature, when only the lowest energy levels are populated, the 
number of 'resonant' levels participating in the absorption increases with increasing 
temperature as T 2 . Since the contribution of each level to the sum l|B.2(l decreases 
as T~ 3 / 2 , the dielectric constant is proportional to i>t- In the opposite case, for high 
temperature, the main contribution is given by the single upper level with f2(e) = u>. 
Therefore, the rate of energy absorption is determined mainly by the derivative of the 
distribution function in the vicinity of this level, which decreases like T~ 3 / 2 ~ v^ 3 
with increasing temperature. 

Collisionless damping may influence reflection and absorption of a laser pulse 
by a film. For a thin film the reflection coefficient is small, R ~ (a/A) 2 <C 1, and 
depends weakly on the imaginary part of the dielectric constant. Therefore, the effect 
of Landau damping on reflection is tiny. We will not consider it here. In contrast, the 
absorption coefficient is proportional to the small parameter a/X and depends strongly 
on the imaginary part of the dielectric constant. To estimate the relative contribution 
of Landau damping to absorption we consider a thin dielectric film irradiated by an 
infrared or optical femtosecond laser with intensity ~ 10 14 W/cm 2 . Laser pulses with 
such parameters are widely using in studies of the optical breakdown in dielectrics 
(see [S] and references therein). The laser field produces inner ionization from the 
valence to the conduction band, where electrons behave as a classical plasma with 
temperature close to the ponderomotive energy in the field 0. 

As follows from equation (|15|l , the energy absorbed per electron in the conduction 
band due to Landau damping is 

87m 



where r is the duration of the laser pulse. For T e ~ 10 eV, a ~ 10 nm and r ~ 100 
fs, from equation (|B.2|I we have x w 0.1 and £l ~ 0.1. The electron density in 
the conduction band usually does not exceed the critical value corresponding to the 
plasma resonance in a dielectric, when u! p — wy'e^ + l. For 1 + e' ~ 2... 3 (i . 
this corresponds to n c ~ 10 21 cm -3 . When the density approaches this critical value 
breakdown sets in. Close to the critical density we obtain from equation (|22|l ol — 100 
eV. This value corresponds to an effective relaxation time for Landau damping of 
tl ~ 10 fs. This is two orders of magnitude less than the electron-phonon relaxation 
time, which is known to be around 1 ps for such parameters. Another mechanism of 

* It should be noticed, however, that such linear behavior is not universal. For example, in the 
'triangular' ID well, U(x) ~ jxj, the dielectric constant is exponentially small at low temperatures, 
so that the law Q is never realized. The dependence of the dielectric constant on the shape of the 
self-consistent potential will be investigated elsewhere. 
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energy absorption is inner ionization. It contributes about 10 eV per electron at the 
same intensity. However, the energy transfer per electron due to Landau damping 
decreases linearly with increasing film thickness a, while the efficiency of absorption 
due to inner ionization practically does not depend on the thickness, at least so long 
as it smaller than the skin depth. Therefore, for films with thickness around 100 nm 
the contributions from Landau damping and inner ionization to the energy absorption 
will be comparable. Obviously, with increasing thickness the effect of Landau damping 
vanishes. 

The main conclusion of the estimations presented above is that for thin films with 
a thickness below a few hundred nm Landau damping may dominate over the other 
mechanisms of energy absorption. As a result, the mean electron energy observed in 
spectra from laser-irradiated films, which is approximately equal to the mean energy 
H22JI absorbed from the field, should decrease proportionally to a -1 with increasing 
thickness a. Indeed, a mean electron energy decrease with increasing thickness was 
clearly observed in experiments with thin carbon films 

The expressions obtained above are valid, strictly speaking, for a monocrystalline 
film. The consequences of Landau damping in polycrystalline films may be different if 
the electrons remain trapped inside their respective crystal layers under the action of 
the laser field. This depends on the relation between the amplitude xq of the electron 
oscillations and the size d of the crystals. If xq <C d then different crystals may be 
considered as independent systems. In this case, the expression for the rate must be 
averaged over the orientations of the crystals so that it will become independent of 
the laser polarization. Besides, spatial averaging will replace the width a of the film 
by the average size d of the crystals. Obviously, this limit is realized for low laser 
intensity (depending on d). In the opposite case of high intensities, when xq ^> d, a 
polycrystalline film is expected to give the same response as a monocrystalline film. 
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Appendix A. Derivation of expression (II 1[) for the absorption rate 

In the (xy) plane, the electron motion is free. Therefore, transitions between different 
energy levels only occur due to the action of the z component of the electric field @ . 
The perturbation operator has the form 

t 

H int (z,t) = — Pz A z {t), A z (t) =-c [ £ z {t)At. (A.l) 
m e c J 

o 

To first order in the laser field, we find 
ire 2 £ 2 

w± = 2 ft m 2^ 2 K n ilPzNI 2 ^ e "i - e » ± M- ( A - 2 ) 

According to the selection rules for dipole transitions, the matrix elements of the 
operator (|A.lf) are nonzero for levels with opposite parity, so that 

m-n = 2k + l, k = 0,±1,±2, ... 
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In the semiclassical limit, we may replace the quantum numbers n and n\ by the 
corresponding energies using the Bohr-Sommerfeld rule and the relations 

? = *7TTT' e ni -e n = nn(e n )(2k+l), p(e) = Rfi(e)/(e), (A.3) 
de Ml{e) 

which can be easily derived from the former. Here f2(e) = 27r/T(e) is the oscillation 
frequency of an electron with energy e in the potential U(z) and /(e) is the 
classical distribution function, normalized to J de/(e) = 1. Taking into account that 
w + (n, ni) — io_(m, n), we may rewrite (|10|l in the form 

Q = HUJ W n,m\P( e n) ~ ( A -4) 

n,ni 

Replacing the sum over the initial state with quantum number n by an integral 
according to 

de 



and performing this integration with the help of the 6 function in equation (|A.2jl , then 
assuming that the typical electron energy is much higher than the photon energy hu>, 
we obtain from equations HA.20 - jA.4jl 

~ Ttn c ae 2 £ 2 z \ - \(ni\p z \n{m))\ 2 d 

Q = —2^T \ m< nM de^^' ^ 

where ra(ni) and e(ni) are functions of n\ owing to the S function in the transition 
probability ljA.2jl . Finally, we introduce a new summation index k instead of n\ 
according to the second relation in (|A.3jl [the resonance condition Ijl2|)] and replace 
the matrix element by its semiclassical limit. For ID systems, the semiclassical limit 
of the matrix element is just the Fourier component of the corresponding classical 
quantity |2"7|: 

T(e) 

(n + s \p z \n)^ Ps (e) = ^ J p z (e, i)e isn ^dt (A.6) 

o 

where p z (e,t) is the classical time-dependent momentum of the electron with the 
energy e. Substituting equation (|A.6j) into (|A.5I) we end up with the rate Ijllll of 
energy absorption. 



Appendix B. The infinitely deep rectangular well with specular and 
diffuse boundary conditions 

Appendix B.l. Specular reflection 

For an infinitely deep one-dimensional rectangular well with width a, we have from 
equations (I5jl. and <|A.6J) : 

. 7r / 2e m c uo 2 a 2 

"(e) =-\—, = 



a\ m e 2n 2 (2k + l 



2 ' 



|p s (e)| 2 = ^f, Z (T)=aJ^, (B.l) 
ir z s z V 27T 
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so that equation (|17fl gives (recall that the superscript s refers to specular reflection) 

, (u^V 128 ^ e- 4 /* 3 ^ 1 ) 2 * 2 / 8T _ i; T 

£l \lu) ^ 6 x 3 f^ Q (2fc + l) 5 ' X V Trm e LU 2 a 2 ~ Loa A ' j 

For a relatively cold plasma in a wide well, we have x <C 1, and equation (|B.2|I 
reproduces the law (JIJ with ?;t in place of «f and C = 2ui 2 /uj 2 . In the opposite limit 
of a narrow well and a hot plasma, x ^> 1, and 

e£«4.13(^V4. (B.3) 

\ cj / x 6 

Let us reminde that the expression l|B.2|) as well as the rate (|A.5|) are valid in classical 
limit, when the condition T 3> tku is satisfied. It is remarkable that for a deep 
rectangular well this restriction can be avoided if (x <C 1). In this case the dielectric 
constant may be obtained directly from the rate l|A.4|) . Using the exact expression for 
the matrix element 

(„ lb » = -4iJ^L I , (B.4) 

replacing both sums by integrals over energies (this is possible only if x <C 1; in the 
general case only one sum can be replaced by an integral) and taking into account 
that e(n) = /2m e a 2 , we obtain instead of (B.2) 

u \ 2 huj 

) xsmhrjKx^), 77 = — . (B.5) 
uj J 21 



Here Ki(ji) is the McDonald function. The areas of applicability of (|B.5J| and (|B.2|I 
overlap: in the classical limit when i] — > 0, the product sinh i]Ki (rj) — > 1 and (|B.5|I 
gives the low-temperature limit of (|B.2ll . 

Beyond the linear approximation, the sinusoidally driven particle in the infinitely 
deep well can serve as a paradigm of chaotic motion The breakdown of the 

invariant Kolmogorov-Arnold-Moser tori for increasing laser field strength has been 
invoked to explain laser- induced damage thresholds 30 . 

Appendix B.2. Diffuse scattering 

In a rectangular well Eq. (|20() with the field I© has the solution 

v+(t, a) = vo (sm(wr + a) — sin a), 

m c uj 

eE . N , , eE 
m c uj 

where the subscript "+" indicates that the trajectory is calculated for an electron that 
starts its motion at the left side of the well at t = (i.e. x+(0) = 0) with velocity 
v(0) = vq = ^j2€/m c when the field has the phase a. If the electron starts its motion 
from the right side, X-(0) — a, then the trajectory is 

V-(t,a) = —Vo (sinful + a) — sm a) , 

m c uj 

/ \ / . . eEp . . . . . , 

X-(t) — a + (—vq H sina i H -(coslujt + a) — cosa). (is. 7) 

m c uj m e uj 2 



x + (t) — (v + — — sina)t+— — -(cos(ujt + a) — cosa), (B.6) 
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The travel times t±{vq,cx) are determined by the conditions x + {t + ) 
x_(r_) = 0. To first order in the field they are 

eE Q 



a and 



a 



it 



cos(z + a) 



cos a + z sin a 



m e u> a 

Then the phase-averaged work Ijl9|l is equal to 



.(B.. 



A(e,r) 



2e e 2 E t 



2 r 



2(1 — cos z) 



(B.9) 



Finally, after averaging over the Maxwellian distribution we get instead of (B.2) 



8ui 



dz 



2(1 — cos z) — z sin z 



-4/ttx 2 z 2 



(B.10) 



where the parameter x is the same as in <|B.2(1 . We can see that in the case of diffuse 
reflection it is not necessary for the electron to satisfy the resonance condition 112JI in 
order to absorb energy: all electrons gain some net energy from the field. However, the 
main contribution into absorption comes from the vicinities of the points z — 7r(2fc+l). 
Using the definition of the variable z in equation l|B.8(l it is easy to find that this 
condition is equivalent to the resonance condition 112fl . 

Finally, this result can be easily generalized to the case when the electron 
experiences exactly one specular reflection before the diffuse one. Physically, this 
corresponds to a thin film with different boundary conditions on different surfaces, 
say to a film having a sublayer. Using the same approach we get 



£t = 



dz 



1 



cos z — z sin z 



(l-cos^e - */* 1 



(B.ll) 



Both for diffuse and for specular-diffuse scattering, in the limit of low 
temperatures (i « 1) the dielectric constant £l goes to zero linearly with x (see 
Fig. 2), regardless of the respective boundary condition. In the opposite limit of high 
temperatures (x ^> 1), it decreases, too, but diffuse reflection leads to a substantially 
lower decrease, proportional to vZ 1 , in comparison with specular scattering. 
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Figure captions 



FIG.l The geometry of the problem. The vector Eo is the electric field of the 
incident laser wave. It is decomposed into its s- and p-polarized components, so that 
E = E s + E p , where the vector E s is parallel to the y axis, while the vector E p lies 
in the plane (nk), where k is the wave vector of the incident laser field and n the unit 
normal vector perpendicular to the plane of film. We use the notation s = E s /E , 
p = E p /E , and $ is the angle between k and n. 



FIG. 2 The imaginary part of the dielectric constant £l as a function of the 
dimensionless parameter x — Vt/ujcl calculated at u = uj p for the ID rectangular 
well with specular (solid line), diffuse (dashed line), and specular-diffuse (dotted line) 
boundary conditions, respectively. The behavior of the dielectric constant at low 
temperatures (x -C 1) is shown in the inset. 



